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SUMMARY 

In this paper a simple two-dimensional model for calculating the rotor- 
dynamic effects of the impeller force in centrifugal compressors and pumps is 
presented. It is based on potential flow theory with singularities. Equivalent 
stiffness and damping coefficients are calculated for a machine with a vaneless 
volute formed as a logarithmic spiral. It appears that for certain operating 
conditions, the impeller force has a destabllzlng effect on the rotor. The 
order of magnitude of this effect can be determined from the stiffness and damp- 
ing coefficients calculated. The paper is a brief review of the author's thesis 
(ref. 14), where more details of the calculation can be found. 

SYMBOLS 

A induced velocity in normal direction in point i from singularity in 

point j 

B Busemann factor 

damping coefficient 
dimensionless damping coefficient 
impeller width 
diffuser width at inlet 
diffuser width at outlet 

absolute velocity of fluid at impeller outlet 
induced velocity 
impeller diameter 
eccentricity of rotor center 
impeller force 

dimensionless impeller force 
stiffness coefficient 
dimensionless stiffness coefficient 
number of elements 
n normal vector 

Q total volume flow in impeller 
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optimtiin volume flow 

total flow relative to optimum flow 
impeller radius 
smallest radius of spiral 
tip velocity of impeller 
velocity from vortex source 

angle between flow velocity and peripheral direction 

blade angle of impeller 

total circulation of impeller flow 

vortex strength per unit length 

density of fluid 


INTRODUCTION 

The linearized governing equation for self-sustained lateral vibration of 
a rotor can be written 

^ + £3. + Kg^ = ^ (1) 

in the absence of external forces. The vector 5 . represents the generalized 
coordinates of the system, ^ is the mass matrix, C the damping matrix, and K 
the stiffness matrix. The stability of the system is determined by the solu- 
tion to the equation 

det(m^ + CX + K) = 0 (2) 

where X = a + i3 and _g^ = Ae^^. If a > 0, the system is unstable- Apart 
from elastic forces and inertia forces due to the deflection of the rotor it- 
self, these matrices also depend on forces from bearings and seals. Further- 
more in turbomachinery, forces from the working fluid acting on the rotor may 
affect the matrices and consequently alter the stability of the rotor (ref. 10 ). 

In centrifugal pumps and compressors the working fluid exerts a force on 
the rotor caused by diffuser/ impeller interaction, as shown experimentally and 
theoretically by different authors (refs. 1 to 9). As indicated by Domm and 
Hergt (ref. 2) and Hergt and Krieger (ref. 8 ), this force depends on the eccen- 
tricity of the rotor. In the literature this force is often called the radial 
force, but this term is a bit misleading in a rotordynamic sense since the force 
also has a tangential component. Therefore in this paper it will be called the 
impeller force . 

In this paper a method is presented for calculating the impeller force, its 
dependency on the rotor eccentricity, and the force's associated stiffness co- 
efficients. Furthermore it is shown that the impeller force also depends on 
the velocity of the rotor center. This gives rise to equivalent damping co- 
efficients, which are also calculated. For simplicity the analysis is re- 
stricted to centrifugal pumps with a vaneless volute but can be extended to any 
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kind of diffuser, as long as neither losses associated with friction and three- 
dimensional flow nor compressibility seriously affects the pressure- 
distribution around the impeller. The calculated stiffness and damping coeffi- 
cients represent the contribution to the stiffness and damping matrices K 
and C from the impeller force and thus permit one to determine the influence 
of the impeller force on rotor stability and synchronous response. 

PHYSICAL MODEL 

This model of a centrifugal pump or compressor stage with a vaneless volute 
is founded on the concept of Csanady (ref. 1). It is based on a two-dimensional 
representation of the diffuser and a representation of the impeller by an equiv- 
alent vortex source concentrated in a single point, as shown in figure 1. 



The flow is considered a potential flow field with the singularity Q,r. 

The vortex-source point coincides with the rotor center, and it may have 
any eccentricity and velocity in a coordinate system fixed to the diffuser. 

The source strength and vortex strength Q and F are related to parameters 
of the impeller as shown later. 

CALCULATION OF FLOW FIELD 

The flow field is calculated by a singularity method, that is, a method 
based on the replacement of the diffuser contour by a distribution of vortices 
and/or sources. The method used was developed by Hess (ref. 11). The diffuser 
contour is considered as a series of small linear segments. Each segment is 
covered by a vortex distribution of uniform strength. The vortex strength var- 
ies from segment to segment. The normal component of the flow field on the 
contour must be zero. If we have N linear segments, we have N vortex 
strengths, and hence we can satisfy the condition in N points. These points 
are chosen as the N midpoints of the segments. 

In a potential flow field with potential singularities the law of super- 
position is valid; that is, the velocity in any point is the sum of the veloci- 
ties induced by all singularities in the field. Thus the kinematic conditions 
give 
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( 3 ) 


N 

Z A_y. = -V .n. 

j=l 

where Vooi is the velocity at the midpoint at section i produced by the im- 
peller flow (through the vortex source F,Q respresentation) with the impeller 
assigned some given eccentricity and velocity. The coefficients are 

evaluated as the normal velocity components induced in point i by a unit vor- 
tex strength along line segment j. The term yj is the unknown vortex 
strength per mit length of line segment j; Vooi is the velocity induced by 
the vortex source F,Q representing the impeller for a given eccentricity and 
velocity of the vortex source. 

The solution of equation (3) gives the values of yj for all j = 1,N. 
After solution of equation (3), all the singularities in the flow field are 
determined, and thus the velocity in any point in the plane can be calculated. 
The development of the coefficients Aj^j and XooilLi shown in appendix A 

for any rotor-center eccentricity and velocity. The method presented is valid 
for a vaneless volute with the rotor center coinciding with the spiral center 
for zero eccentricity. It could be developed for any diffuser geometry as long 
as the pressure distribution around the impeller is satisfactorily calculated 
by means of potential theory. Equation (3) will normally contain a large num- 
ber of linear equations - from 100 to several hundreds - and the coefficients 
of the system are generally different from zero. 

DETERMINATION OF Q AND F FROM MACHINE DATA 

Figure 2 illustrates the velocity vectors of the flow at the exit of the 
impeller. With Csanady’s nomenclature (ref. 1) the variables of this figure 
are defined as follows: 

B Busemann slip factor, tabulated by Wislicenus (ref. 12) for different 
impellers; depends on blade angle and number of blades 

B blade angle of impeller 

U tip velocity of rotor 

C absolute velocity of fluid at outlet 

The total volume flow is 


Q = 2Trr^bj^ sin a^C (A) 

where r^ and b^ are the radius and the width of the impeller, respectively. 
From figure 2 

C(cos Uq + sin aQ)cot B = UB (5) 

Combining equations (A) and (5) yields 

Q = UB27Trj^bj^/ (cot + cot B) (6) 
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Furthermore 


Q = tan a^-r ( 7 ) 

because Oq is the angle between the flow and the peripheral direction at the 
outlet. The optimum efficiency with a spiral-formed volute with spiral angle 
X is obtained when 


at the volume flow Qopt> according to Csanady (ref. 1). From equation (6) we 
obtain 


tan a^(tan 6 tan X + 1) 

^rel “ ^^^opt " (tan 3 + tan a^) 

CALCULATION OF IMPELLER FORCE 

The force per unit width on a body with circulation T in the presence of 
a source Q in a parallel stream with the velocity c is, according to the 

theorems of Joukowski and Lagally (in ref. 13) ~ 

^ = p (jc Q + ^r) (9) 

where p is the density of the fluid. If we consider the velocity induced by 
the singularity distribution on the diffuser contour to be a parallel stream in 

the calculation of the impeller force, the Impeller force is obtained by in- 

serting equations (6) and (7) into equation (9) to yield 

UB2iTr^bj^ 

-i = 1 + tan OQ/tan S “o^i^ 
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In this equation, is the velocity induced in the rotor center by the con- 

tour singularity distribution, as determined from equation (3). The solution 
procedure for ^ for any rotor-center eccentricity and velocity is carried 
out in appendix B for a vaneless volute formed as a logarithmic spiral, whose 
center coincides with the rotor center for zero eccentricity. The procedure 
requires the vortex distribution solution yj from equation (3). This could 
be done for any diffuser geometry. 

CALCULATION OF EQUIVALENT STIFFNESS 
AND DAMPING COEFFICIENTS 

The impeller force can be calculated for any eccentricity and velocity of 
the rotor center as shown on the previous pages. The force is now calculated 
for different combinations of eccentricity and velocity of the rotor center in 
the X- and y-directions of figure 3. 



Figure 3 


The values of the force can then be tabulated or stored in the computer , and 
the equivalent stiffness and damping coefficients are calculated as the follow- 
ing numerical derivations: 


K. . = -AF./Ae 
ij i 


B. . 

ij 


-AF./Av. 
1 J 


(i=l,2;j=l,2) 


( 11 ) 


where the index 1 corresponds to the x-direction, and the index 2 corresponds 
to the y-direction. 


RESULTS 

On the following pages some calculation results are shown for a given im- 
peller with different volute-spiral angles. The force and stiffness and damp- 
ing coefficients are presented in the following forms: 
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DISCUSSION 


From figures 5 and 6 we see that for certain relative flows the cross- 
coupling stiffness coefficients KXY and KYX have opposite signs. This means 
that the impeller force has a destabilizing effect on the rotor. Furthermore 
the damping coefficients BXX and BYY become negative for certain values of the 
relative flows (figs. 11 and 12). This negative damping is highly destabiliz- 
ing. For the cross-coupling damping coefficients BXY and BYX, we have almost 
exactly BXY = -BYX (figs. 9 and 10), and these coefficients are 4 to 10 times 
larger than BXX and BYY. This means that the rotor is subject to relatively 
large gyroscopic forces. These forces would tend to stabilize the rotor. 

The spiral angle X is a significant design parameter, as seen in fig- 
ures 4 to 12. For increasing X the impeller force decreases in magnitude, 
the stiffness coefficients increase in magnitude, BXX and BYY decrease in magni- 
tude, and BYX and BXY increase in magnitude. The interval where the destabiliz- 
ing effects are absent is moved to the right on the Qrel ”^xis for increasing 
X as seen from figures 5, 6, 9, and 10. The effects of impeller design parame- 
ters are not investigated in this paper. This effect could, if wished, be cal- 
culated by equation (8). 

The only relevant measurements reported in the literature are those deter- 
mining the impeller force as a function of the relative flow (refs. 2 to 5 and 
9) . There is qualitative agreement between the reported values of the magni- 
tude of the impeller force and the present calculations. As for the direction 
of the force there is a large scatter in the results reported in the literature, 
and the results of this paper consequently only agree with some of these re- 
sults. No direct measurements of stiffness and damping coefficients are re- 
ported in the literature. Experimental determination of stiffness and damping 
coefficients associated with the impeller force must be made before a compari- 
sion of calculated and measured values is possible. 
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APPENDIX A 


CALCULATION OF INFLUENCE COEFFICIENTS A.. AND V n FOR A 

ij — »i— i 

VANELESS VOLUTE FORMED AS LOGARITHMIC SPIRAL 

For this purpose we have to use some basic formulas from potential-flow 
theory. The velocity induced in a point (x,y) by an infinately long line 
source in (0,0) is 


(v)q = l/(x^ + (13 

where Q is the source strength. For a vortex line of strength T, the corre 
sponding velocity is 

(“) - + y2)(_y) (14 

If we have a line segment placed as in figure 13, covered with a constant vor- 
tex strength of unity, an integration over the segment gives the following 
velocity induced by the line segment in point (x,y) (ref, 11): 



V 

X 



V 

y 


“In 


(x -I- As/2)^ 
(x - As/2)^ 



1/ Att 




J 


(15) 
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where Vy ^ 0 and Vx ->■ 1/2 for x = 0, y ->■ 0. Consider the logarithmic spiral 
of figure 14, 



Figure 14 


Take 6^ and as the values of the angle 6 for the position vector 

of the midpoint of the line segment i and j. We have from figure 14 


n , = 


cos(0^ + X) 
-sin(0^ + X) 


= r^e j 


„ , , /sin 0 . 

0 . /tan X/ j 


cos 0 


sin(0j + X) cos(0j + Xyi/Xj 


-cos(0^ + X)-sin(0j + 


/x^\ r-sin(0. + X)-cos(0. + X)|/x - ^in 0 A 

lyW cos(0^ + X)-sin(0j + X) W - ^cos 0^^ 

These coordinate transformations are necessary for the calculation of the 
induced velocities from the vortex distribution. Furthermore we have for the 
induced velocities in the different coordinate systems 

/''xlnX . (19! 


-cos(0j + X)-sin(0j + ^y’ .in 


If the point (x^.y^^) is the midpoint of line segment i, we have 


= r^e J 


0./tan i 


cos 0 . 
>. 1 
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Hence, by using the coordinate transformation (18), the coordinates of point 
(Xi»yi) in reference system (xj,yj) are 



_e0i/tan ^ ~ ^cos X 

i/tan “ ®i^ ^sln X 


( 21 ) 


From equations (16), (21) and (15), we get , the normal velocity Induced in 

point i from line segment j 


Ay -0 


for i = j 


Ij 4 tt 


-sln(0^ - 6j)2 


X. . + As/2' 
arctan I — ; 

y. . 

ij 


'x:. - As/2^ 
- arctan \ — — i 

y . . 

ij 


+ cos(0 - 0 .)ln 

i J 


( x^^ + As/2)^ + y^^ ^ 
,(^ij - As/2)^ + y_ 



for 1 j 


( 22 ) 


In addition, from 


n.V . 


-“1— 


equations (13), 

“cos(0 . + X) 
1 

sin(0^ + X) 


(14), (7), and (20) 
^ Q/2irx^^(r^ 


, we obtain 

- cot ot^f.) 
0^1 


with 



This is the induced velocity in the normal direction in point i from the 
vortex source (Q,F) in the impeller center with eccentricity (e^,6y) and veloc- 
ity (v^jVy). For this geometry a correction is required because the volute 
does not have a constant width. Domm and Hergt (ref. 2) use the correction 
factor t)Q/b2, where bQ is the width at inlet and b2 the width at outlet of 
the volute, so the final result is 


n.V . 


— 1—001 



-cos (0 , 

+ X) 


/ V \ 

*'0/»2 

1 

sin(0 . 

+ X) 


1 1+ Q/2Trr^^ (^^ - cot oi^^) 


L 1 



\y/ 


(23) 
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APPENDIX B 


CALCULATION OF INDUCED VELOCITY IN ROTOR CENTER 


From equations (15), (18), and (19), we can calculate the induced velocity 
contribution for each line segment in any point in the plane. If we calculate 
these contributions in the rotor center, the total induced velocity in the 
rotor center is the following sum from all N line segments of the volute con- 


tour: 



But the rotor center has the velocity (v^.Vy) , so the induced velocity is 



(24) 
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